typeset using JPSJ.sty <ver.l.Ob> 



(N 
O 
O 
(N 

> 
O 



Role of Collective Mode for Optical Conductivity and Reflectivity 
in Quarter-Filled Spin-Density- Wave State 

Yuh TOMioB and Yoshikazu SuzumuraQ 

Department of Physics, Nagoya University, Nagoya 464-8602 
(Received ) 

Taking account of a collective mode relevant to charge fluctuation, the optical conductivity of 
spin-density- wave state has been examined for an extended Hubbard model with one-dimensional 
quarter-filled band. We find that, within the random phase approximation, the conductivity ex- 
hibits several peaks at the frequency corresponding to the excitation energy of the commensurate 
collective mode. When charge ordering appears with increasing inter-site repulsive interactions, 
the main peak with the lowest frequency is reduced and the effective mass of the mode is en- 
hanced indicating the suppression of the effect of the collective mode by charge ordering. It 
is also shown that the reflectivity becomes large in a wide range of frequency due to the huge 
dielectric constant induced by the collective mode. 

KEYWORDS: optical conductivity, collective mode, spin-density-wave, charge ordering, reflectivity, quarter-filled 
band, extended Hubbard model 
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§1. Introduction 

It has been well known that coUective modes of the 
charge fluctuation play an important role to determine 
the property of density waves, e.g., the pinning of the 
density wave and the dynamics^ siiidi. as the optical con- 
ductivity in organic conductor.S'clEroEl' 

The effect of the collective mode of density waves on 
the optical conductivity is found in the pronounced peak 
at zero frequency, which gives rise to a sliding motion 
of the density wave. The collective mode of the charge- 
density-wave (CDW) state leads to the zero-frequency 
peak of the conductivity, whose magnitude depends on 
the strength of the coupling to phonon. For the spin- 
density-wave (SDW) state, the conductivity is deter- 
mined only by the collective mode where the full weight is 
given by the peak. Such a fact is well known only for the 
incommensurate density wave. The optical conductivity 
for incommensurate density waves has been extensively 
examined since the work by Lee, Rice and AndersonEP 
who calculated the conductivity by taking account of the 
collective mode for both CDW and SDW. The collective 
mode, which participates in the sliding motion of the 
CDW, has the effective mass enhanced by the electron- 
phonon interaction. For large effective mass, the dom- 
inant contribution to the conductivity comes from the 
single particle excitation. On the other hand, the con- 
ductivity of SDW state exhibits only the zero frequency 
peak with the unrenormalized electron mass, i.e., the 
contribution from the single particle excitation dimin- 
ishes due to complete compensation with that of the 
collective mode and the collective mode determines the 
conductivity for all the frequencies. The conductijirity of 
SDW has been calculated analytically by TakadaE^ The 
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numerical calculation of the conductivity with the sev- 
eral electron-phonnn interactions has beep performed by 
Fenton and Ares.B Virosztek and MakicP examined the 
pinning effect and the optical conductivity by introduc- 
ing the phase Hamiltonian. 

The calculation of the conductivity for the commen- 
surate case is complicated since all the higher harmonics 
of the density wave must be taken into account. For un- 
derstanding the conductivity in the presence of the com- 
mensurability and interaction, it is useful to calculate a 
sum rule showing that the summation of the conductiv- 
ity with respect to the, ftequency is equal to the average 
of the kinetic energyEj'U^ The sum rule depends on the 
magnitude of the interaction except for the limit of the 
weak coupling of the interaction as seen from the incom- 
mensurate case. The exact numerical calculation for the 
extended Hubbard model at quarter-filling has shown the 
reduction of the averaged value. E3' This indicates an en- 
hancement of effective mass of density wave. The optical 
conductivity of the Hubbard model is further examined 
by introducing the dimeriaation, which has an effect of 
producing the charge gap.tS However the calculation of 
the optical conductivity based on the mean field theory 
is required for the SDW state since the exact calculation 
of one-dimensional system leads to the absence of the 
long range order. 

The SDW state of organic conductors has been studied 
extensively in terms of the mean-field theory. In order to 
understand the experiment_ij|_Drganic conductoES-where 
SDW coexists with CDW,yEJ several authorfO'B-EJ 
have examined the ground state of the extended Hub- 
bard model at quarter-filling. They have shown such a 
coexistence when the inter-site repulsive interaction be- 
comes larger than the critical value. Based on the mean 
field ground state, the collective modes for charge fluctu- 
ation have been calculated where the charge gap of the 
collective mode vanishes at the onset of charge order- 
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This comes from the second order phase tran- DmQo = ^ ^ (^Cl^Ck+ 



mg 

sition between two kinds of ground states. Thus it is of 
interest to examine a role of such a low lying excitation 
on the dynamical quantities. The optical conductivity 
at quarter-filling was calculated within the single parti- 
cle excitation,LiP and the effects of collective excitations 
on absorption spectra were studied by applying the ran- 
dorrL^jikase approximation to the systems with a finite 
size.cHIE3' The purpose of the present paper is to demon- 
strate such a role of collective mode for the optical con- 
ductivity and the reflectivity. Further we examine the 
effect of charge ordering, which may originate in intcr- 
site repulsive interaction. 

In §2, formulation for the optical conductivity is given. 
Based on the mean-field ground state, the random phase 
approximation is applied to calculate the response func- 
tion for the electric current. In §3, the optical conduc- 
tivity is calculated numerically. The effective mass is 
estimated from the main peak of the conductivity. In §4, 
the effective mass of the collective mode is calculated in 
terms of the phase Hamiltonian, which is derived from 
the response function. The result is compared with that 
of §3. Based on the optical conductivity, we examine the 
dielectric constant and refiectivity in §5. The plasma 
frequency is calculated from the sum rule of the opti- 
cal conductivity. The results is analyzed in terms of the 
collective mode with low lying excitation. In §6, the ef- 
fect of dimerization on the conductivity and the effective 
mass is briefly discussed. 

§2. Formulation 

We consider a onc-dimcnsional extended Hubbard 
model given hy H ~ Hq + Hmt, where Hq and iJint are 
respectively defined by 



N 



h.c. 



N 

Hint = E (^"-iT^ii + VrijUj+i 
i=i 



V2njnj+2) 



(2.1a) 
(2.1b) 



In the Hamiltonian, Cj^. denotes a creation operator 
of an electron at the j-th site with spin a = (1,1), 
Tij — rij^ + Tiji and rijc — Cj^Cja- A periodic bound- 



ary condition is taken, i.e., C| 



Cj^ with the to- 



tal number of lattice site N. The quantity t denotes 
the transfer energy and corresponds to the dimeriza- 
tion where t and the lattice constant are taken as unity. 
Quantities U, V and V2 are the coupling constants for 
repulsive interactions of the on-site, the nearest-neighbor 
site and the next-nearest-ncighbor site. 

For a quarter-filled band, where the Fermi wave num- 
ber is given by kp — 7r/4, the mean-fields (MFs) repre- 
senting SDW and CDW are written as (m = 0, 1, 2 and 
3, and sgn(o-) = +(-) for a =t (i)) 



-E 



E 

T,i -7r</c<Tr 



Sgn(cr) (C^^Cfc+^Qo 



MF 



(2.2a) 



(T = t,X —TI<k<7T 



MF 



(2.2b) 



where Qo{= 2fcF = 7r/2) induces density waves with a 
periodicity of four lattices. The z-axis is taken as the 
quantized axis of the SDW ordered state. The expres- 
sion (0)mf ^fipotes an average of O taken by the MF 
Hamiltonian,!!^ which is expressed as 



MF 



NU 



J2 5]*Ufc)i?MF(fc)*.(fc) 



0<k<Q„ a 
\---- 

V 16 2 



Qo\ 



■-\Sc 



- - D 



■'■2Q0 



-S. 



2Qo 



+ NV{-=:+DIq^ 
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+ NV2{-- + 2m 



2Qo 



(2.3) 

with 5'J-(fc) being the four component vector given by 

Hik) -/C^L,C't^Qo,.,Cl+2Qo,.'^ri+3Qo,J- TheMF 
Hamiltonian has the following matrix elements. 



Hupik) 



-21 cos(fc + mQo) + U/A + V + V2 



\ / m+l,m Z Z 



Mf(^) 



m,m+l 



(2.4a) 



(2.4b) 



(HuFik)) = i^-2V+2V2)D2Q,-sgn{a)^S2Q, 

\ / m-t-2,m Z Z 

+ 2iid sin(fc + mQo) . (2.4c) 



In the [Present case of quarter-filling, MFs can be rewrit- 
ten aM> So = 0, - 1/2, Sq, - S*^q^^ = Sie'', 
Dqo = ^3Q„ = D,e^^o--f'\ S2Q0 = 5*Q„ EE S2 and 
D2Q0 = D*Q^ EE D2. Quantities Si{> 0), L»i(> 0), 
52 and D2 denote amplitudes for 2fcF SDW, 2A:f CDW, 
4fcF SDW and 4A:f CDW, respectively. The quantity 
6 denotes a phase of SDW. The phase diagram of the 
ground state on the plane of y_and. 14 has been shown 
explicitly in the previous paper OEj'cH' There are three 
kinds of ground state on the plane of V and V2 where re- 
spective order parameters are given by S'l in the region 
(I), S'l and D2 in the region (II) and Si, S2 and Di in 
the region (III). The region (II) (region (HI)) is given for 
large V {V2). 

The MF Hamiltonian, eq. (2.3), is diagonalized with 
the eigenvalues Ei{k) (Z=l,2,3 and 4) and eigenvector 
\lka) where four energy bands are defined as Ei(k) < 
E2{k) < Ej,{k) < Ei{k). These quantities are given 

by, 



FMF(fc) \lka) ^ Ei{k)\lka) 

\lk(7)„i = Fi„ra[k) , 



(2.5a) 
(2.5b) 



where Fima{k) is the component of the eigenvector. We 
note that only Ei{k) is occupied due to quarter-filling. 
The effect of commensurability of quarter-filling is suffi- 
ciently large in the present choice of J7 = 4, where the 
gap between Ei{k) and E2{k) is much larger than the 
band width of Ei{k). 
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(a) 



(b) 




Fig. 1. Feynman diagram for the conductivity ,0) where each lines 
consist of the basis of four components. The solid and wavy lines 
denote the electron Green function and the effective interaction, 
respectively and the cross denotes the current vertex. Respective 
diagrams denote (a) current-current response function for single 
particle excitation [ito('»^)], (b) current-current response function 
for collective mode [crc(aj)] and (c) the effective interaction. 



In order to examine the optical conductivity, we use 
the following lattice version.ll2P The current operator is 
expressed as 

J = ieJ2{t- i-iyt^) - Cj,Q+i., 

JO" 

3 



0<k<Qo,<J m=0 



ium{k) Cl 



+mQo,cr^k+{m+2)Qa,a 



(2.6) 



where Vm{k) — 2tsin(A:-|-m(5o) ^^nd u,n{k) — 2fdC0s(fc-|- 
mQo). In terms of eqs. (2.6), the Kubo formula for the 
dynamical conductivity is given by 



N 



aiu;)+ia'{u;) = - / dt e'"* ([<z(t), J(0)]) 



(2.7) 



where <t{uj) and cr'(w) are real and q = ej^ja J ^ja-^jo-- 
The optical conductivity is given by cr{uj). Within the 
random phase approximation, the real part of eq. (2.7) 
is calculated aaj^ 



cr(cj) = cro(w) + ac{uj) 



(2.8) 



where cro('^) a-nd CTc('^) are respectively given by Fig. 1(a) 
and 1(b) as explained below. For the convenience of 
numerical calculation, we introduce a damping factor rj 
by a> — > w + iyy. Thus a delta function corresponding to 
a pole is replaced by a peak with a finite width. 

The first term of eq. ( ^.8[ ) denotes the contribu- 
tion from single particle excitation, which is written as 
[Fig. 1(a)] 

ao(c.)=Re||-[n°^.(^)-nO/0)]| , (2.9) 

where 

U%ic.)^^iyT{TrJir)jm 



f{Ei,ik))-f{Ei{k)) 

0<k<Qo, ami, 1712 1,1' ' ^ ' ' ' 

X [Vmiik)Fi„i^cr{k) - iUrm{k)Fi^„n+2,a{k)] 



E EE 



X [vinJk)Fi*^^„{k)+m„^,{k)Fi*^^^2,Jk)] 

X Fi.^2.{k)Ft,„,^Ak) , (2.10) 

and Tt is the ordering operator for the imaginary time 
T. The quantity /3~^ denotes a temperature with the 
Boltzmann constant taken as unity and ujn — 2Trnf3~^ 
with integer n. The function f{x) denotes a Fermi dis- 
tribution function given by [e(="-^)^ + l]-i where the 
chemical potential is given by ^ = {Ei{Qo) +i?2(Qo))/2 
in the limit of zero temperature. In deriving eq. (2.10), 
we used a single particle Green function (the solid line in 

Fig. 1), written as (G"^(fc, iwi.))"^ = iuj^- HMpik) +M 
where is the Matsubar a fre quency for Fermion. 

The second term of eq. (j^^ denotes the contribution 
from the collective mode expressed as 



a,(c^) =Re<j -n^=^.(a;) 



(2.11) 



where 



= ^r t {TrJ{T)Jm e'""-| 



2N 



niu(w)- 



U{0) 



1- u{0) n(o,w) 



-n,y(w) . (2.12) 



The matrix, JJ (g)/(l— U (q) U {q,iLj)), denotes an ef- 
fective interaction (Fig. 1(c)). The quantity, [/((?), has 
only the diagonal element given by [U {q))ii' = 5ii'Ui{q), 
{i = 1,3,0,2) where Ui{q) = UM - Uo{q) - 
-(L/-l-4t/ cos 9-1-4^2 cos 2g) andi72(g) = -{U-'iVcosq+ 
4V2COs2g). In eq. (2.12), respective vector and matrix 
are written as 

n,„(c.) = (n^"i(c.),n°f(a;),n°S(^),n^"^(^)) (3.13) 
n„,(c.) = (ni°(a.),n3;^(a.),n°^(^),n|(.;))*, (2.14) 



( nll{q,uo) ni3(g,a; 

nS^(g,-) n°3 
V n2i(g,^) nl(g 



n3o 



23/ 



nf,('?,^) 
n°/,(9,c.) 
nf,(9,c.) J 



(2.15) 

The elements of eqs. (2.13)-(2.15) are calculated as 



n°r(c.') -^jyr{TrJ{r)SUO)) 



MF 



I iuJn—*uj-\-iri 



-A. V V V f{Ei,{k))^f{Ei{k)) 
2N ^ ^ LO+iri+Eiik)-Ei{k) 

0<k<Qo,ami,m2l,l' ' ' ' 

X [Vmi{k)Fi.mia{k) - iUmi{k)Fi^„n+2,a{k)] 

X Fi',m+ni2,<7ik)Fi'm2aik)Fi*,^^„{k) sgn(cr) , 

(2.16) 



— dr(T.5.„.(r,0)J)Mpe'--|^^^_^j^ 

V V- V- f{Ei,{k))-f{Ei{k)) 
2N ^ ^ to + 171+ El, (k)~Ei(k) 

0<k<Qo,<jmi,m2 1,1' ' ^ ' ^ ' 
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X Vr. 



X Fi 



(2.17) 



and 

r/3 



E EE 



/(i;p(fc))-/(i;;(fc + g)) 

27V ^ ^ uj + iri+Ei'(k)-Ei{k + q) 

a<k<Qo,<ymum2 Ll' ' 'W U 

XF, {k + q)F*„^,^„^^^^ik + q) 

X Fi'^^„{k)Fi%^^Jk) sgn(cr) , 



(2.18) 

where ^^^(g) = Efe^r ^L^*:+9+mQo,CTSgn(cr) and 
Dmiq) = Y^ka^L^k+q+mQo.a- Quantities ifjTi'^)^ 
n2J°(w) and n^™'(g,a;)(= n™™'((j,w)) are given by 
eqs. (2.16), (2.17) and (2.18) without sgn{a), respec- 
tively. 

The spectrum of the cohective mode is obtained from 
det(l- Uiq) U{q,io{q))) -0 , (2.19) 
where the basis vector is given by0-' 

Hq) ^ {S,iiq),SM'l),Mq),D2{q)y . (2.20) 

The sohition of eq. ( ^.19[ ) consists of several poles, i.e., 
isolated solution. We define ujc{q) as the pole with the 
lowest energy which corresponds to the mode leading to 
the sliding motion of density wave io. the limit of weak 
interaction (or incommensurate casecP). In addition to 
Wc(0) expressing the charge fluctuation, some of other 
poles with higher energy may also give rise to the peak 
of the optical conductivity if they describe the charge 
fluctuation. In the next section, these mode are exam- 
ined by rewriting the basis of eq. ( 2.2C ). In eq. ( 2.19 ), 
there are also three continua due to quarter-filled band, 
which contribute to the conductivity. 

For calculating the effec tive mass, we note sum rules 
for a{uj) and adt^) of eq. (2.8) expressed as 



TO r- 

= — %/2( 



L 



dui crc(w) = , 



(2.21) 



(2.22) 



where to(= 7r/(4\/2)) is the band mass, JLn these equa- 
tions, the former is known analyticallyllJ'til' and the lat- 
ter is verified numerically in the present paper. The 
mass, TOc, which is reduced to m for the incommensu- 
rate case, is renormalized by interactions. The effective 
mass of the collective mode can be obtained from the 
weight of the main peak of the conductivity. 

§3. Optical Conductivity 

The frequency dependence of optical conductivity is 
examined to comprehend the role of the collective mode 
of the commensurate density wave at quarter-filling. Pa- 
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Fig. 2. The density of states for U = 4, = and with the fixed 
(a) V = V2 = (region (I)), (b) V = 1.6 and V2 = (region 
(II)), and (c) V2 = V/1.2 = 1.4 (region (III)). 



rameters U, V and V2 are chosen as intermediate cou- 
pling e.g, U — 3, 4: and 5. We take = in this section 
and discuss the effect of 7^ (i.e., dimerization) in the 
final section. 

Since the fundamental contribution to the optical con- 
ductivity comes from the single particle excitation, we 
calculate the density of states defined by 



Im{TrG"'(fc,cj + iO-^)} , (3.1) 



^ ^ ttN ^ 

0<fc<Qo 

i.e., D{cj) = 2/{nN)\ims^+o Ek,. Ei '^/((^ " Ei{k) f + 
5^) where / dujD{Lj) ~ 2. In Fig. 2, the density of states, 
D{u!) is shown for U = 4, where three choices of parame- 
ters of V and V2 correspond to region (I), region (II) and 
region (III), respectively. The occupied states are shown 
by the shadow, the area of which is 1/2. At quarter fill- 
ing, density of states are divided into four parts, corre- 
sponding to Ei{k), E2{k), E^{k) and Ei{k). There exists 
a symmetry with respect to w = (i?2(Qo) + E3{Qo))/2 
only for F = V2 = 0. The gap between the first band 
{Ei{k)) and the second one (£'2(fc)) increases with in- 
creasing V and/or V2 while the second gap is nonmono- 
tonical due to vanishing a,t V = U /A in the region (II) as 
seen from eqs. (2.4b) and (2.4c) with 6 = {]. 

Now we calculate the optical conductivity. Figure 3 
shows the conductivity for the Hubbard model with U 
= 3, 4 and 5, where dotted curve and solid curve cor- 
respond to oqIuj) and cr(w), respectively. The dotted 
curve shows the three continua, which correspond to the 
excitation from the filled (i.e., first) band to the sec- 
ond band, the third one and the fourth one. The solid 
curve exhibits three sharp peaks, which have the follow- 
ing relation to the spectrum of collective modes aX q — Q. 
For [/=4, the peaks are located at ~ 0.4, 3.0 and 4.8 
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Fig. 3. Optical conductivity (7{lo) is shown by the solid curve for 
U = 3, 4 and 5 with V = V2 = 0, = and r] = 0.02 where 
(j{w) = (Jo{u!) + ac{u>) and the dotted curve denotes (tq, the 
contribution from single particle excitation. 



Fig. 4. Optical conductivity u{lu) is shown by the solid curve for 
V = 0, 0.4 and 1.6 with (/ = 4, V2 = 0, td = and r] = 0.02 
where the notations are the same as Fig. 3. 



while eq. ( 2.19 ) with q ~ gives five poles at uj — 0.38, 
2.24, 3.05, 3.88 and 4.80 respectively with three continua 
0.91 <uj < 2.00, 2.42 < uj < 3.01 and 3.91 < to < 4.42. 
It is found that the peaks of (t{uj) correspond to the first, 
third and fifth poles. The contribution from the second 
and fourth poles is absent within the numerical accu- 
racy of the present calculation. For U = 3, the second 
and third peaks of cr{Lu) are reduced. This indicates a 
fact that <t{llj) in the limit of small U is reduced to that 
of the incommensurate case where all the contribution of 
cr(uj) is given by the first peak with cj — > and the contri- 
bution from the single particle (i.e. aolcu)) compensates 
completely with that of CTc(w). With increasing U, the 
weight for the lowest peak decreases and the second and 
the third peaks increase together with the continuum es- 
pecially in the frequency region just below the second 
peak. 

The effect of V on cr(w) is examined in Fig. 4 for U =4 
and V2 = where <t{uj) with V = 0.4 and 1.6 are com- 
pared with that oi V — 0. The conductivity of the up- 
per panel is the same as that of ?7 = 4 in Fig. 3. The 
region (II) is obtained for V > Vc{= 0.34) while the 
region (I) is obtained for V < Vc. For V — 0.4 (the 
middle panel), the peaks are located at w ~ 0.2 and 
4.0 while eq. ( |2.19| ) with 9 = gives four poles at lu 
= 0.18, 2.21, 3.98 and 4.65 respectively with three con- 
tinua in the regions 0.92 < cu < 2.11, 2.28 < w < 2.93 
and 4.00 < uj < 4.48. Thus, the peaks of a{uj) corre- 
spond to the first and third poles. The effect of CTc is 
small for intermediate frequencies. This may come from 
the fact that the gap between the second band and the 
third band (i.e., £^2(0) and £^3(0)) vanishes a,t V = U/A 
in the region (II) as seen from eq. (2.4c) and 6 — 0. For 
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Fig. 5. V dependence of m/rric (dashed curve) and mc/m* (solid 
curve) for U = 4, V2 = and = 0. The ground state in the 
region (I) [region (II)] is obtained for F < Vc(~ 0.34) ( V > Vc). 



V = 1.6 (the lower panel), the peaks are located at a; ~ 
1.7, 2.9 and 4.7 while eq. ( |2.19| ) with q^O gives five poles 
at UJ = 1.67, 1.70, 2.89, 2.96 and 4.71 respectively with 
three continua 1.79 < uj < 2.54, 2.97 < uj < 3.56 and 
4.74 < a; < 5.04. Thus, the peaks of a{uj) correspond to 
first, third and fifth poles. It is found that, within the 
numerical accuracy, the second and fourth poles of the 
collective mode do not give the peak of the conductivity. 
With increasing V, the weight of the first peak is re- 
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CO 



Fig. 6. Optical conductivity cr{Lo) is shown by the solid curve for 
V2 = 0.8 and 1.4 with U = i, V = I.2V2, td = and r; = 0.02 
where the notations are the same as Fig. 3. 




duced. For large V, the compensation of the continuum 
by (Tc(w) is reduced and the continuum of (7(0;) becomes 
noticeable. 

Here we examine the total weight of the cond uctiv ity. 
In Fig. 5, the quantity m/rric defined by eq. ( 2.21 ) as 
the function of V is shown by the dashed curve where 
TTic increases with increasing V due to narrowing of the 
filled band. Note that m/nic = 1 in the limit of small U 
and V. Here we examine the effect of V on the weight 
of the main (i.e., first) peak of cr(a;). Since the finite 
imaginary part rj is introduced by treating w — > w + ijy, 
we calculate the weight using a method given by 



Wo 

du! fTc(l^) 



(3.2) 



where Uq denotes the lowest value of uj satisfying aduj) = 
0, i.e., the solid curve and the dotted curve in Fig. 4 cross 
each other at oj = loq. The quantity m* may be regarded 
as the effective mass for the collective mode with the 
first (lowest) pole. In Fig. 5, the quantity rric/m*, which 
denotes a fraction of the collective mode is shown by the 
solid curve. With increasing V, the quantity mc/ni* de- 
creases indicating the decrease of the contribution from 
collective modes. There is a jump of mc/m* at V ^ Vc- 
In the limit of the weak coupling, one obtains rric/m* = 1 
showing that the conductivity is determined only by a 
collective mode with w ^ 0. 

The effect of V2 on the conductivity is examined in 
Fig. 6 with [/ = 4 and V/V2 = 1.2. Note that the param- 
eter with V2 — 0.8 (V2 = 1.4) corresponds to the region 
(I) [region (HI)] as seen from a fact that the boundary 
between region (I) and region (II) for U — A is given 
byllf V V2+ 0.34 for V2 < V2c{^ 1.32). Since the 
ground state of V2 — 0.8(= V/1.2) is the same as that of 
V — V2 — 0, the dotted curve of V2 = 0.8 is the same as 
that of F = in Fig. 4 (and also [/ = 4 in Fig. 3). How- 
ever the solid curve deviates from that oi V = due to 
fTc(w), which depends on not only U but also V and V2. 
The peaks correspond to the first, third and fifth poles 
of the collective mode where their five poles located at 



Fig. 7. V2 dependence of m/rric and mc/m* for U = 4, V = 
I.2V2 and t(j = 0. The ground state in the region (I) [region 
(III)] is obtained for V2 < V2ci^ 1.32) ( V > V2c)- 



UJ = 0.29, 2.24, 3.03, 3.88 and 4.72 are slightly different 
from those of Fig. 3 with U — 4. Such a small change 
of a{uj) is due to the effect of V2 on the charge fiuctua- 
tion. For V2 = 1.4, there appears five peak s located at 
u ~ 0.9, 2.7, 3.2, 4.1 and 4.8 while Eq. ( |2 19)) with q = 
gives six poles at = 0.85, 2.59, 2.71, 3.25, 4.12 and 
4.80 respectively with three continua 1.15 < ui < 2.01, 
2.72 <uj < 3.23 and 4.15 < a; < 4.59. Thus, the peaks 
of (t{uj) correspond to the first, third, forth, fifth and 
sixth poles. We make following remarks as the effects 
of large V2. The collective mode except for the second 
pole gives rise to the peak in the conductivity. The con- 
tinuum of (j{uj) which exists in the intermediate range of 
frequency is mainly determined by cto(w), i.e., that of the 
single particle excitation. The contribution from the sin- 
gle particle excitation becomes dominant with increasing 
V2. 

Using the result of Fig. 6, we examine the V2 de- 
pendence of rric/m* and m/mc, which are defined by 
eqs. (2.21) and (|3.2|). In Fig. 7, rric/m* and m/nic as the 



function of V2 are shown by the solid curve and dashed 
curve respectively for U = A and V/V2 = 1.2. It is rea- 
sonable that m/rric is constant for V < V2C, i.e., region 
(I) since the quantity given by eq. (2.21) is determined 
only by the ground state. The quantity m/mc (dashed 
curve) exhibits a cusp a,i V — V2C but that of Fig. 5 
is continuous at V — Vc- This fact is comprehended 
as follows. The ground state of region (HI) is obtained 
by a competition between spin density (5'i) and charge 
density (Di) while that of the region (II) (i.e., the coex- 
istence with D2) is obtained mainly by the change of the 
phase of SDW. It is noticed that rric/ m* is independent of 
V2 for V2 < V2C although the pole of the collective mode 
depends on V2. The rapid increase of m* for V2 > V2C 
shows a fact that the charge ordering {Di ^ 0) induced 
by V2 enhances strongly the effective mass of the collec- 
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tive mode. For the case of varying V (the sohd curve of 
Fig. 5), such an enhancement takes place for V, which is 
much larger than Vc ■ Thus it is found that charge order- 
ing of Z?2 is less effective compared with that of Di. 

Here we note the relation between the peak of the con- 
ductivity and the pole of the collective mode. As men- 
tioned in the previous section, the collective mode gives 
rise to the peak of the conductivity if it is optically ac- 
tive, i.e., the mode can describe the charge fluctuation. 
It is well known that, in the mode of the incommensu- 
rate case, the components of 2fcF SDW and that of — 2fcF 
SDW) are the same in the magnitude but opposite in the 
signfl-' Extending the sta te to the present case of quarter- 
filling, we examined eq. (|2.19|) in the limit of small q by 
choosing a restricted two-component base as 

^q) ^ (e-^'S^q) - e^'SMq), D2{q)f . (3.3) 

We have verified that all the pole obtained from eq. (3.3) 
correspond to the location of the peak of the conductivity 
in Figs. 4 and 6. For the commensurate case, there are 
many poles describing the charge fluctuation while the 
pole with lowest energy gives the dominant weight. 

§4. Effective Mass of the Collective Mode 

The effective mass to* of the collective mode has been 
estimated from the weight of the main peak of the con- 
ductivity, in §3. In this section, we examine the 
effective mass of the collective mode with lowest excita- 
tion energy using a different method and compare with 
the result of the previous section. For this purpose, we 
utilize the phase. HamiltoniarLJiihich has been studied 
for both COWsEI) and SDWs.BN) Using phase variable 
9{x) which describes the charge fluctuation of SDW with 
long wave length, the effective Hamiltonian is written as 



n^(^) 



de{x) 

dx 



bqt,e\x) 



(4.1) 



where [9{x),I{p{x')] — \5{x — x'). In stead of deriving 



parameters, pflj b and q^ in eq. (4.1) from path inte- 
gral method,EZP we determine them numerically using the 
charge susceptibility and the spectrum of the collective 
mode [eq. ( 2.19D ], which has a gap due to commensura- 
bility. 

The coefficient, 6, is estimated as follows. From 



eq. (4.1), the Green function is calculated as 

_ 1/2 



/ dT(T,0,(T 

Jo 



)^_,(0))e" 



(4.2) 



where 9q is the Fourier transform of 9{x). It is noted that 
the Fourier transform of charge fluctuation with wave 



number q is given by iq9q/'K. Then, from eq. (4.2), the 
static charge susceptibility for small q, is written as 

Xc(g)-(J)'(V-,> = i-^ . (4.3) 



27r26 g2 



9o 



The charge susceptibility is also calculated from the re- 
sponse function, which is deflned bycil^ 



2N 



At (T,$(g,r)$t(g,0))e" 



+a;-|-iO 



(4.4) 




Fig. 8. The static charge susceptibility Xc as a function of q/Qo, 
for V = 0, Vc{~ 0.34) and 1.6 where U = A, V2 = and = 
0. In the inset, the numerical results of x<^{l) (dots) for the 
small q are compared with solid curves obtained from Xc(9) = 
g2[27r2b(g2-h<7g)]-l. 



where $(<?) is given by eq. (2.20) and the response func- 
tion is evaluated in the random phase approximation as 



- n(g,w)/(l- U{q) Ii{q,oj))- From eq. (^, the static 
charge susceptibility is also obtained as 

xT^{q) = {Do{q)Do{-q)) , (4.5) 

where Do(q) = iq9q/Tr. Thus the quantity b can be 
obtained by comparing eq. ( [4.3| ) and eq. (4^). In Fig. 8, 
Xc{q) as the function of q/Qo is shown for V — 0, Vc 
(dotted curve) and 1.6. One finds Xc(0) — except for 
V — Vc which leads to Xc{q) 7^ for q ^ 0. The results 
ior V = Vc is reasonable if we note eq. (^) does hold 
and the collective mode becomes gapless, i.e., qo = 0. In 
the inset, Xciq) obtained from eq. ( [4.3[ ) (solid curve) is 
compared with that of eq. (4.5) (dots) where b in eq. (4.3) 
has been determined to fit the dots. The quantity a is 
calculated by substituting b into the excitation spectrum 



obtained from eq. (2.19), which can be written as i-iJc{q) — 
sjbiq^ -\- qQ)/a for small q. In the following calculation, 
we use renormalized quantities defined as a = Anvpa and 
b = ^-nb/vp where vp — ^/2. It is found that a — b — 
1 in the limit of weak interaction. We note that d = 



m*/m since the first term of eq. (4^) may lie rewritten 
as TTVF{m/m*)Il'j, in a way similar to CDWO Thus the 
quantity d can be compared with the effective mass of 
the previous section, i.e., eq. (3.2). 

In Fig. 9, quantities d and b as the function of U 
are shown by the solid curve. The open circle is ob- 
tained from eq. (3_^) where ij = 0.02 for J7 > 4 and 
rj = 0.02{E2iQa/2) - Ei{Qo/2)) otherwise for numer- 
ical accuracy. The open circle coincides well with d. 
There is the visible enhancement of the effective mass 
m* for 11^4:. The U dependence of b is also shown 
by the solid curve. For comparison with 6, we show 
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Fig. 9. (7 dependence of a(= 47ri)pa) and b{= 47rfe/?;p) for V = 
= = n wh ere a and 6 are coefficients of the effective Hamil- 
tonian eq. . The dashed line denotes 1 + IfJJnvp ) . The open 
circles represent m* jm estimated from eq. ( p.^ ). In the inset, 
the solid curve and the doifeted curve are Kp{= 1/ (ab)^/"^) and 
the exact one, respectively.E3' 
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Fig. 10. V dependence of d (solid curve) for (7 = 4 and V9 = Ia = 
0. The open circles represent m* /m obtained from eq. (3.2). In 
the inset, solid curve and closed circle denote Kp and the exact 
result 



1 + U/{ttvy) by the dashed hne, which is the formula 
in the Umit of smaU U and V = V2 = 0- The soHd 
curve deviates from the dashed Hne for large U due to 
the effect of commensurability. The good coincidence of 
h and the dashed line for small U indicates that eq. (4.1) 
describes well the effective Hamiltonian as seen in one- 
dimensional Hubbard model.tj' In the inset, we show 
Kp{= l/(a6)^/^), a parameter for the charge fluctuation, 



Fig. 11. V2 dependence of a for C/ = 4, V = I.2V2 and = 0. 
The open circles represent m* /m obtained from eq. (3.2). 



by the solid curve where the dashed curvej-denotes exact 
one for one-dimensional Hubbard model.EB' The result 
indicates a fact that the charge fluctuation around the 
mean-filed ground state may be justified even quantita- 
tively. 

Further we examine V and V2 dependence of a. In 
Fig. 10, the normalized effective mass, m* /m, as the 
function of V is shown by the solid curv e wh ere the cir- 
cle denotes the result obtained by eq. (3.2). Both re- 
sults remain almost the same and independent of V for 
y ^ 1. The mass enhancement is found for V{^1) which 
is much larger than Vc- The enhancement of m* obtained 
from the response function is smaller than that obtained 
from the optical conductivity. In the inset, V dependence 
of Kp is shown by the solid curve where the closed circle 
denotes the exact result H It is found that the charge 
fiuctuation obtained by the present approximation may 
be justified also in the presence of V . 

In Fig. 11, V2 dependence of the normalized quantity 
m* /m is shown where open circle denotes the result ob- 
tained by eq. (3.2). Both are almost the same and stay 
constant for V2 < V2c- In contrast to Fig. 10, m* /m in- 
creases rapidly for V2 > V2c- The effective mass obtained 
from eq. ( ^.2[ ) is larger than that of the present section. 
The former is the direct calculation while the latter is 
the result under the assumption of phase Hamiltonian of 



eq. (4.1). The reason for the discrepancy for large V and 
V2 remains unclear at present. 

§5. Dielectric Constant and Reflectivity 

In order to examine the reflectivity, we calculate the 
dielectric constant e(w), which is given by 

e{oj) = l + — {a{uj)+ia'{uo)) , (5.1) 

UJ 

where <j{uj) + icr'(w) is the dynamical conductivity de- 
fined by eq. ([2.7]). The imaginary part, (t'(w), can be 
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calculated using the Kramers-Kronig relation with the 
optical conductivity cr{io). However we calculate <j'{uj) 
directly from eqs. ( |2.9| ) and (2.11). Actually, one obtains 
(t'((jj) = (Jq{uj) + (t'c{cj), where 



(5.2) 



(5.3) 



In Fig. 12, 
and V = V2 



cr'((jj)/e^ and Re[e(a;)] are shown for U ~ 4 
~ ~ 0. The quantity a'{Lo)/e^ is the odd 



function with respect to uj. One finds that a'{uj)/e oc 
—Lo for small w and changes the sign at u; = Wc(0). With 
increasing w, Re[e(a;)] increases rapidly but decreases to 
change the sign at the frequency, i^c(O), corresponding 
to the main peak of (j{oj). The lower panel at a; = 
reads Re[eo(0)] = 1 (dotted curve) and Re[e(0)] ~ 62 
(solid curve). Thus, the effect of the collective mode 
is noticeable even at low frequency. This fact can be 
comprehended by the phase Hamiltonian derived, in §4. 
From eq. (4.2), the conductivity is calculated 

ia'{uj) 



cr(w) 



(5.4) 



where 'D{q,iuj. 
r.h.s. of eq, 
the dielectric constant is estimated as 



is the Green function defined by the 
( f4.2| ). Substituting eq. ( [s!^ ) into eq. (|]l]), 



1 - 



(5.5) 



where Wc(0)(= qa^b/a) corresponds to the pole of the 
collective mode with lowest energy. The quantity u-p 
denotes a plasma frequency defined by 



It is found that eq. 
Re[e(ti;)]. We note that 



explains 
e(^) = 1 - 



(5.6) 

the behavior of 
{uj'^rac/m)/uj^ in 



the limit of large u where ujprric 

In ter ms of the complex dielectric constant given by 
eq. (|3.l|) , the reflectivity is calculated as 



1 



(5.7) 



In Fig. 13, R{uj) is shown for U ^ 4, V2 = U = 
and fixed V — 0, 0.4 and 1.6 where the solid curve is 
obtained from the conductivity with the collective mode 
and the dotted curve is obtained from the conductivity 
without the collective mode. The arrow denotes the loca- 



tion of the plasma frequency defined by eq. (5.6) where 
we take e = 1. The upper panel corresponds to the 
result for the conventional Hubbard model with the in- 
termediate interaction. The dotted curve exhibits the 
insulating behavior (i.e., small R{i-u)) not only below the 
gap but also in a certain region above the gap where the 
cusp of the dotted curve at a; ~ 0.9 corresponds to the 
gap for the single particle excitation. The solid curve 
shows that the collective mode enhances the reflectivity 
in a wide region above the gap of the collective mode 



3^ 



O' 



o- 

o- 
o 

<N - 







CO 



Fig. 12. u) dependence of the imaginary part of the conductivity 
[o-'(t.ij)] and the real part of the dielectric constant [e(t<j)] for U = 
A, V = V2 = and t;j = where the solid curve (dotted curve) 
has been evaluated from cr' (o"q). 




Fig. 13. ui dependence of reflectivity R{lo) for U = 4, V2 = and 
t(j = with the fixed V = 0, 0.4 and 1.6 where the solid curve 
(dotted curve) has been evaluated from a + ia' (ctq -|- icr' ) . The 
arrow denotes a plasma frequency ob tain ed fro m c o. (5.6). The 
dashed curve is calculated from eqs. ( ^.e| ) and (5.7). 



[wclO)]. This comes from the huge magnitude of the di- 
electric constant as seen from Fig. 12 and eq. (5.5). Here 
we note the following fact indicating a significant role 
of the collective mode with Wc(0). The dashed curve 
in Fig. 13 is calculated from eqs. (|5.5|) and (5.7). One 



finds a good coincidence between the solid curve and the 
dashed curve in the wide frequency region except for high 
frequency regime where the effect of commensurability 
becomes complicated. Thus the reflectivity in the most 
region of frequency is determined by the collective mode 
with the lowest excitation energy. The sharp dip just 
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Fig. 14. uj dependence of reflectivity R(yi) for [/ = 4, 1/ = I.2V2 
and t(j = with the fixed V2 = 0.8 and 1.4 where the notations 
are the same as Fig. 13. 



below the plasma frequency is found for the frequency 
close to the second peak of the conductivity (first panel 
of Fig. 4). The dip at w = 3.03 originates in the fact that 
Re[e(a;)] takes a local maximum (~ 0.13) at a; = 3.03 and 
becomes negative for a; < 3.00 and 3.04 < w. The effect 
of third peak of the conductivity is invisible in the re- 
flectivity. The collective mode enhances the reflectivity 
at low frequency but suppresses it at high frequency. We 
note that the solid curve becomes the same as that of 
the free particle (i.e., the metallic state) in the limit of 
small U , which leads to the conductivity oc l/(iw — r?). 

Now we examine the effect of inter-site repulsive in- 
teraction on R{iS). The middle panel of Fig. 13 shows 
the behavior for V just above where the charge order- 
ing begins to develop. Compared with the upper panel, 
i?(w) is enhanced at low frequency due to the decrease of 
the gap of the collective mode. The small dip just below 
ojp originates in the single particle contribution while the 
sharp dip at w ~ 4 comes from the collective mode. We 
note the special case of F = 14 which leads to the gap- 
less excitation of the collective mode. In this case, i?(w) 
is expected to be similar to that of free particle except 
for the region at high frequency just below the plasma 
frequency as seen in the middle panel. The lower panel 
shows the behavior for V much larger than where the 
charge ordering develops well. We find the insulating be- 
havior of R(iS) in the wide range of low frequency and 
still i?(a;) < 1 above the gap. Thus, i?(w) exhibits the in- 
sulating behavior for all the frequencies when the charge 
ordering well develops to localize the electron. Since the 
difference between the solid curve and the dashed curve 
increases with increasing V , the charge ordering reduces 
the effect the collective mode with cjc(O). In Fig. 14, 
R{ijS) is shown for = 4, VjVi = 1.2 and = 0, and 
fixed V2 = 0.8 and 1.4 where the notations are the same 
as Fig. 13. The upper panel is similar to that of 1/ = 
of Fig. 13 since both ground states are in the region (I). 
The lower panel is the result for the ground state in the 
region (III) where charge ordering comes from V2- The 
insulating behavior is also found but is slightly compli- 



T 1 1 1 1 1 1 1 1 1 r 
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2 4 6 

(0 

Fig. 15. Optical conductivity (y(ij^) is shown by the solid curve 
for t(j = (solid line), 0.08 (dotted line), 0.16 (dashed line) and 
0.24 (dash-dotted line) with (7 = 4, y = V2 = and ?? = 0.02 
where the notations are the same as Fig. 3. 



cated at high frequency. Two kinds of dip are located 
close to the second and third peak of the conductivity 
(lower panel of Fig. 6) where the origin is the same as 
that of Fig. 13. Based on Figs. 13 and 14, it is found 
that the reflectivity is suppressed sufficiently by charge 
ordering. 

§6. Discussion 

We have examined the effect of the collective mode 
on the optical conductivity and the reflectivity for SDW 
with one-dimensional quarter-filled band using the ex- 
tended Hubbard with intermediate coupling constant. 
The collective mode with lowest excitation energy ex- 
hibits a dominant contribution not only for the limit of 
small interactions but also for the intermediate coupling. 
When U increases to the order of band width, in addi- 
tion to the main peak of the conductivity, other peaks 
become visible at high frequencies, which correspond to 
poles of the mode being active to the charge response. 
The effect of the collective mode decreases in the? pres- 
ence of the inter-site interaction, which induces charge 
ordering. The effective mass estimated from the weight 
of the main peak of the conductivity increases when the 
charge ordering is induced by inter-site repulsive interac- 
tions being larger than the critical value. The enhance- 
ment of effective mass is also analyzed in terms of the 
phase Hamiltonian. Both results agree well for weak in- 
teractions while the latter seems to underestimate in case 
of strong interactions. The effect of the collective mode 
on the reflectivity is further examined. We found that 
the mode with lowest excitation energy determines the 
reflectivity in a wide range of frequency between the sin- 
gle particle gap and the plasma frequency. The reflectiv- 
ity in the range remains large for SDW without charge 
ordering while it is reduced in the presence of charge 
ordering. 

Here we examine the effect of dimerization which is 

introduced in cq. (2.1a). The main effect of id on the 
ground state is to extend the region (I) into region (II) 
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Fig. 16. td dependence of m/m* (solid curve), mc/m* (dotted 
curve) and raj vac (dashed curve) for V = 4, V = V2 = and 
r] = 0.02. 

(the region (III) is less effected)^-' since the dimerization 
is compatible with the spatial variation of spin density 
of the region (I). The effect on the charge gap of the 
collective mode is rather complicated. The charge gap 
vanishes at the critical value, corresponding to the 
onset of charge ordering, which increases with increasing 
ti,x3 Thus id enhances the charge gap for small V^(< 
while td suppresses the charge gap for large V . In the 
rigorous treatmen±_.of one-dimensional Hubbard model 
with dimerization,Ej' the dimerization is crucial to obtain 
the main (coherent) peak. In this case, the location of 
the peak is shifted from zero frequency to finite frequency 
since the charge gap is produced only by dimerization. 
In the present case where SDW forms the long range 
order, the charge gap does exist without dimerization 
indicating a quantitative effect of on the coherent part. 
However one may still expect a noticeable effect on the 
continuum (incoherent) part at intermediate, frequency 
as does for the pure one-dimensional case.E3 Actually 
we obtained the following results. Figure 15 shows the 
conductivity with td = 0, 0.08, 0.16 and 0.24 for C/ = 4 
and V = V2 — where the upper panel denotes (Tq and 
the lower one denotes a. With increasing t^, (Tq moves 
to the regime of high frequency. For cr, one finds several 
characteristics. The weight of the first peak decreases as 
the function of id and a shoulder just above the first peak 
appears. The second peak decreases (increases) for td < 
0.16 (0.16 < td) and the third peak decreases (increases) 
for td < 0.08 (0.08 < td) with increasing td. Another 
peak appears just below the second peak. These come 
from the interplay of the excitation of the collective mode 
and that of the single particle. In Fig. 16, the effective 
mass is shown as the function of td. The total weight 
of the conductivity {mc/m) remains almost constant or 
even slight decreasing as the function of td while the 
effective mass {m* /m) determined by the first pole is 



increased by td. 

Finally, we comment on the relevance of the present 
result to the optical experiment in the materials with 
the one-dimensional quarter-filled band. Although the 
experiment exhibits various kinds of additional effects 
such as the coupling to phonon and one-dimensional fiuc- 
tuation in the presence of the correlation gap, it is ex- 
pected that, in organic conductors, the small refiectivity 
of TMTTF salt compared with that of TMTSF salt could 
be due to the presence of charge ordering. EP 
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